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Abstract 

A general multi-step N i — ► M probabilistic optimal universal cloning pro- 
tocol is presented together with the experimental realization of the (1 — * 3) 
and (2 — > 3) machines. Since the present method exploits the bosonic nature 
of the photons, it can be applied to any particle obeying to the Bose statistics. 
On a technological perspective, the present protocol is expected to find ap- 
plications as a novel, multi-qubit symmetrizator device to be used in modern 
quantum information networks. 
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A most relevant limitation in quantum information processing is the impossibility of 
perfectly cloning (copying) any unknown qubit \<f>) [1]. Even if this process is unrealizable in 
its exact form, it can be approximated optimally by the so-called universal optimal quantum 
cloning machine (UOQCM), one which exhibits the minimum possible noise for any possible 
input state. From a theoretical perspective, two different kinds of universal cloning machines 
have been developed so far: a deterministic N — > M UOQCM based on a unitary operator 
acting on N input qubits and 2(M — N) ancilla qubits [2,3] and a probabilistic UOQCM 
based on a symmetrization procedure involving a projective operator acting on N inputs 
and (M — N) blank ancilla qubits [4]. 

In the last years several experimental realizations of the UOQCM for polarization (n—) 
encoded photon qubits have been reported. The deterministic UOQCM has been realized 
by associating the cloning effect with QED stimulated emission [5] while the probabilistic 
machine has been realized adopting a linear symmetrization protocol [6] . Thus far, only the 
simplest 1 — > 2 cloning processes, i.e. for N = 1 and M — 2, were realized by both schemes. 
In particular, the probabilistic process was achieved by exploiting the bosonic character of 
the photons within a linear Hong-Ou-Mandel interferometer scheme [7] . 

The present work presents the first generalization of the universal optimal cloning process 
by the realization of a very general linear procedure to extend the probabilistic protocol to 
any value of N and M according to a suggestion by Werner. [4]. The validity of this 
theoretical scheme is supported by the here reported experimental implementations of the 
1 — > 3 and 2^3 probabilistic processes for 7r— encoded photon qubits (n— qubits). 

Let us outline first the N — > M probabilistic cloning theory. Consider N identically 
prepared unknown qubits in the state pi = \<p) (<p\ as input of the cloning machine while 
(M — N) blank qubits, i.e. all in the state Pa — \i are used as an auxiliary resource. To 
generate M output clones the machine performs the symmetrization of the output state by 
applying the projective operator, 11+^ over the symmetric subspace of M qubits: 
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where Pn^m = 2 m-n j^jf is the success probability of the procedure. All the identical 
output clones are described by the same output density matrix a N ^ M = Tn-^m \<f>) (<t>\ + (1 — 
Fn^m) 1^) (0% where T N ^ M = (<f>\<J N ^ M \<P) = (N+l+P)/(N+2) with (3 = N/M < 1 is 
the "fidelity" of the optimal cloning process [2,3,8]. In the case of our present experiment, 
the value of these parameters for the 1 — > 3 and 2^3 machines are found respectively: 

n 1 77 _ 7 „ _2 r _ 11 

Pl^3 — 21 •'1^-3 — gj £>2^3 — 3, •/ 2->3 — if- 

In order to implement the process expressed by Eq.l, consider that any generic 1 — > 
M cloning process can in facts be realized by a chain of (M — 1) intermediate identical 
machines according to the operatorial identity: n^=n^(IT^ M_1 ' ) <g> I^ 1 ^)=II^ M ' ) (II^ M ~ 1 - ) ® 
l(i))(n^ M_2) ®l( 2 )) • .(ir+ 2) <g> l( M ~ 2 )). This is justified by the very definition of the symmetric 
subspace of M qubits as the smallest subspace in Hf M spanned by the tensor product vectors 
\<p)® M for any \<f>) G H 2 , the qubit space with d=2. In the above expression iW stands for 
the identity operator in the i— qubit space, Hf \ This implies that the symmetrization of 
M qubits can be carried out step by step e.g. starting from the symmetrization of the two 
input qubits p and p, as shown in Fig. la [8]. Precisely, the state realized at the output 
of any i th machine in the chain, i.e. of the overall 1 — > % device, belongs to the set: 

e &> = nf(p ® p)nf, .. e w = nJV*- 1 ) ® p)n« ..^ = nf V (M_1) ® p)nf } = 

n (M ) (n (M-l) _ (n (2) 0n (M-2) )[p0 -«(M-l) ](n (2) 0l (M- 2)) . ^ ® I«)nf } (2) 

Note that in the above expressions the input states, i.e. the pure p = |0) (0| and the 
mixed p = |, can be interchanged leading to the two equivalent configurations shown in 
Fig. la: the upper one has been chosen for the present implementations. These schemes 
are represented by arrays of equal Hong-Ou-Mandel interferometers, each one consisting 
of a 50/50 beam-splitter (BS) and realizing the qubit symmetrization in Hilbert spaces 
of increasing dimensions. The theory above can be easily extended to the analysis of any 
general N — > M cloning process. Optionally, the procedure could be made to consist of a 
sequence of linear symmetrization devices acting by inequal cloning steps, e.g. by injection 
of different mixed states along the chain. 
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The Fig. lb shows the experimental apparatus by which the 2-step chain (1 — > 3) = 
(1 -> 2) + (2 -> 3) UOQCM has been realized. 

1 — > 2 UOQCM. The device realizing the first step state-symmetrization was the beam- 
splitter BSa excited over the two input modes by the 2-qubit state: p},^ 2 = \4>) (<f>\ s ® ^f. 
After projection by BSa in the symmetric subspace, the output state realized on the output 
mode k 2 was expressed by: 

Plut = Q {2) = 1 100) (001 + 1 1{00 ± }> <{00 X }| (3) 

where the notation [{00^}) stands for a total symmetric combination of the states \4>) and 
\4>~ L ). The identical condition realized on mode k\ was neglected, for simplicity. The two 
clones j = 1,2 emitted over k 2 were expressed by the same operators: o- l -^ 2 =Tr h ^jp 1 Q ^ 2 = 

110) (0I + ||0 ± ><0 ± |- 

In the experiment a pair of photons was generated by spontaneous parametric down 
conversion (SPDC) in a 1 mm thick BBO crystal, cut for Type I phase matching. The two 
photons, each with wavelength (wl) A = 795nm and coherence time r co h — 200/s, were 
emitted over two modes fc^ and kg respectively in the product state of horizontal (if) linear 
polarizations (it): \H) S \H) A . Then, on mode ks the qubit \H) S was n— encoded by an 
optical waveplate (wp) WP^ into the generic pure-state \<f>) s , ps = |0) (01^, while on mode 
kA the qubit \H) A was transformed into a fully mixed-state pa = \ by a depolarizing channel 
realized by a stochastically driven Electro-Optics Pockels (EOP) cell, Pa- The two qubits 
Ps and pa were then drawn into a linear superposition in BSa- The exact space-time overlap 
of the two input modes implying the actual realization of the interference was controlled by 
the microscopic BSa displacement: Z A = 2cAt Let's call " BSa- interference" the condition 
Za = corresponding to maximum interference. By turning on the cloning machine, i.e. 
setting it in B S a- interference, the induced Bose coalescence implied an enhancement by a 
factor i?i_».2 = 2 of the \(fxf>) component in the 2-qubit output state and no enhancement of 
the [{00 -1 }) component [6]. The measurement of was carried out by a post-selection 
technique, by the n— analysis setup shown at the r.h.s. of Fig. lb, connected directly to 
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the output mode k 2 , by disregarding at this stage the presence of BSb- This n— analyzer 
consisted of an output mode selector, realized by a 5 meter long single-mode optic fiber, 
followed by the wp WP^ 1 that mapped the output state \<f>) into \H) by counterbalancing the 
action of the input WP^. Finally, by a polarizing-B S ', (PBS) the \H) and \V) components 
of the output state were directed respectively over the modes k$ and The mode was 
coupled to the detectors D x , D 2 , D 3 by means of two equal 50/50 beam-splitters BS 1 , BS 2 
while the mode k^ was coupled to £>*, D 2 by BS 3 . The \<p<p) component was identified by 
detecting coincidence by the D— pair sets [D±, D 2 ] and [D 1 , D 3 ] while the state [{(/x^}) was 
identified by [Di,Dl\. The detectors (D) were equal single-photon counters SPCM-AQR14. 

Three different input states \<j)) s = \H), 2-*(\H) + \V)), 2~*(\H)+i \V)), identified in the 
following by \H) , \H + V) and \H + iV) respectively, were adopted to test the universality 
of the device. The cloning process was found to affect only the \(fxf>) component, as expected, 
and R\^ 2 was determined as the ratio between the peak value (cloning machine switched 
on) and the basis value (off). The corresponding experimental values of the cloning fidelity, 
F^ 2 = (2.Ri^ 2 + 1) / (2.R1-2 + 2) were: F£ 2 = 0.831 ± 0.001; F?+? = 0.833 ± 0.002; 
F\^ 2 W = 0.830 ± 0.002. These values are to be compared with the theoretical value F*^ 2 = 
5/6 0.833 corresponding to the optimal enhancement ratio R = 2. Similar results for the 
1 — > 2 UOQCM have been reported in [6]. 

1 — > 3 UOQCM. In agreement with the upper configuration shown in Fig. la, the 50/50 
beam-splitter BS B was the next state-symmetrization device: the whole (1 — > 3) UOQCM 
is shown in Fig. lb. This BS was excited over the input mode k 2 by the output state 
Pout 2 °f tli e i — * 2 UOQCM, Eq. 3, and over the other input ks by the fully mixed state 
Pb = ^f- In analogy with the first step experiment, this state was obtained by means 
of a EOP, Pb acting on a highly attenuated quasi single-photon beam expressed by the 
n— qubit \H) B , deflected from the main laser by the mirror M, and delayed by Z B = 2cAt B 
via an "optical trombone". Once again, the condition Z B = 0, dubbed here as "B SB- 
interference" condition, was made to correspond to the maximum overlapping of the 2 
input modes of BSb- In particular, in no-" B S 'a- interference" condition, i.e. for \Za\ 3> 
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2cr co h, the " B SB-interference" corresponded to the maximum overlapping in BSb of the 
mixed states pa and p B . In summary, the overall 1 — > 3 machine was excited by the 
input state p}^ 3 = |0) {4>\ s <S> ^<8> \, i.e. by the pure state ps to be cloned and by two 
mutually uncorrelated mixed states pa and p B . By applying to this state the projector 

nf = (000| + ({000^}) <{000 X }| + |{00 ± ± }> <{00 ± ± }| + |0 ± ± ± > <0 ± ± ± |, 

the symmetrized output state is obtained: 

p^ 3 = = jj |000) (000| + | |{000 ± }> <{000 X }| + 1 |{00 ± ± }> <{00 ± ± }| (4) 

Each one of the identical clones j = 1,2,3 can be thought to be expressed by the reduced den- 
sity matrix: aj^ 3 =Tr h)k ^jpl~^f=l |0) (0| + | |0^) (0^1- The projection over the symmetric 
subspace was identified by the measurement of the 3-photon Fock state over the output 
mode k 3 by the 7r— analyzer apparatus already described and shown at the r.h.s. of Fig. lb. 
The output field emitted over k 4 was negleced. The |000), [{000^}) and |{00" L 0" L }) com- 
ponents of Po^ 3 were measured by the 3-fold coincidence events respectively by the detector 
(D) sets [D U D 2 ,D 3 ], [D u D m ,D* n \ and [£>,, D* n , D*} for any l,m = 1,2,3 and n,p = 1,2. 
A little inspection of the circuit leads to the following expectations. In the exact B SB- 
resonance, and no-B S A-resonance, i.e. in the condition "Bose coalescence" of only the two 
mixed states p^and p B in BSb, an enhancement by a factor T of the |000) and |{00~ L ± }) 
components should be detected by the n— analyzer. Furthermore, by turning on also the 
B S ' A-resonance, i.e. by setting Za= Zb= 0, a further enhancement of the |000) compo- 
nent by a factor i?^ 3 and of the |{000 ± }) component by a factor R\_+ z were expected. In 
summary, the full resonance condition, corresponding to the swithing on of the pro- 
jector, implied the global enhancements by the factors TR\^ 3 , and V respectively of 
the components |000), |{000~ L }) and |{00 ± 0~ L }) of the state p^ 3 , Eq.4. Accordingly, the 
first step of our strategy consisted of the measurement of Y. This was provided by injecting 
in the apparatus the pure state \^) SAB = I^Os \H) a W)bi by setting Za= Z b = and by 
turning off the mixing EOP devices Pa and P B . The value of T was determined by the 
ratio of the counting rates of the |{000 ± }) components, the only non vanishing one under 
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\^)sab excitation. These rates were measured at the peak of the detected resonance curves, 
i.e. with Z S b = 0, and far from the peak, with Z SB » 2cr coh , being Z SB = Za — Z B 
the mutual delay between qubits S and B at BS B . The measured value T exp =1.66±0.05, 
expressing the degree of indistinguishability attained between photons coming form different 
sources, SPDC and attenuated laser, was to be compared with the theoretical one V th = 2. 
By restoring the full operation of the overall apparatus under excitiation by p}^ 3 , R\_> 3 
and -Rf^ 3 were determined as the ratios of the \<j><f>4>) and component measured, 

via 3-D coincidences, in resonance, i.e. Z B = Za = and out of resonance, i.e. Z B = 0, 
Za» 2cT co h- From these measurements, the fidelity of the overall process could be deter- 
mined: F±^ 3 = (3I\RJ_ 3 + 4R?_ 3 + r)/(3r J R^ 3 + 6i2f_ 3 + 3r). The plots of Fig. 2a show 
the experimental 3-D coincidence results measured by the n— analyzer for the various state 
components of the output state: p^ 3 . The experimental values of the fidelity measured 
by the above procedure under excitation of three different input states ps = \(f>) {(f>\ s were: 
= 0.758 ± 0.008; F^ 3 V = 0.761 ± 0.003; FfL% v = 0.758 ± 0.008 to be compared with 
the optimum value: F^ 3 = 7/9 « 0.778. 

The protocol can be easy generalized for any long UOQCM chain by a straightforward 
repetition of the above procedure, as follows. The output state p^ 3 of the 1 — > 3 UOQCM 
be injected into one input arm of a further state-symmetrizing 50/50 beam-splitter BSc 
while a mixed one-photon state pc = be injected on the other arm. This state is generated, 
as previously, by extracting by a further mirror M a highly attenuated beam expressed by 
the 7T— qubit \H) C and mixing it by an EOP, Pc- This will result in a 1 — > 4 UOQCM 
apparatus generating the output state p^ 4 .Then again: the output state p^ 4 of the 1 — > 4 
UOQCM be injected into one input arm of a state-symmetrizing BS B while a mixed state 

p D = 2^. ...and so on. 

2^3 UOQCM . As a significant variant of the above protocol, the Pa Pockels Cell was 
removed and the beam-splitter and BSa was excited over the input modes &u and ks by the 
same pure states: qa = Qs = 10) (<t>\- In BSa— resonance condition, i.e. Za = 0, the BSa 
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acted as a conventional Hong-Ou-Mandel interferometer emitting over the output mode fc 2 
the symmetric Bose state \<p) (0|f 2 to be injected, together with p B = |l B into the BS B 
according to the discussion above: Fig. lb. The input state to this novel N — > M UOQCM 
with N = 2 and M = 3 was then expressible as: pf^ 3 = |0) (0|f 2 <8> ^f, to be mapped onto 
the 3 clone output state, according to: 

Pin 3 - Pl7t 3 = \ 1000) (0001 + \ |{000 X }> ({000^} | (5) 

Each output clone j = 1,2,3 could be expressed by the state o-j^ 3 =Tr h:kf Ljpl^ t 3 = y| \<j>) (<f>\ + 
i^|0 ± ><0 ± |- 

Once again, the \4>4>4>) and [{000^}) components of p 2 ^ 3 were measured by the 3-D coin- 
cidence events respectively by the sets [D 1: D 2 , D 3 ] and [D h D m , D*] for any I, m=l,2,3 and 
n = 1,2. The plots shown in Fig. 2b express the experimental results. In full analogy with 
the previous discussion, the detected bosonic coalescence was expressed by an enhancement 
of the detected coincidences by a factor -R 2 _+ 3 = 3 of the \4>4>4>) component of pf^ 3 whereas 
no enhancement affected Furthermore, in analogy with the 1 — > 2 cloning pro- 

cess, the experimental value of i?2^3 was determined as the ratio between the peak resonant 
value for the BS B interferometer (Z B = 0) and the out of resonance value: Z B » 2cr. The 
universality condition was assessed by injecting the same three different input test states 
adopted for the previuos cases. The experimental values of fidelity were found : F^ 3 = 
0.895 ± 0.003; F^ 3 V = 0.893 ± 0.003; F^ 3 iV = 0.894 ± 0.003 in correspondence with the 
input states defined above. Tese values are be compared with the calculated optimal value: 
F*t 3 = (3i? 2 ^3 + 2) / (3i? 2 ^3 + 3)= (11/12) « 0.917. In all previous experiment, the mea- 
sured values of i?2^3 were reduced by the unwanted injection of two and three photons in 
the mode k B and by simultaneous emissions of two pairs from SPDC. These spurious events 
affected the measured value of the i?2^3 factor by a calculated average amount of 15%. 

In summary a very general and efficient linear multi-step optical procedure for the prob- 
abilistic N — > M optimal universal cloning machine has been proposed together with the 
successful experimental realization of of the first two steps, i.e. the (1 — > 3) UOQCM. Fur- 
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thermore, the probabilistic (2 — > 3) UOQCM, was also demonstrated by a straightforward 
variant of the same protocol. This shows that a very similar protocol can be adopted to 
implement contextually the N -> M UOQCM, the N — > (M - jV) Universal NOT gate 
and any programmable anti-unitary map [9] following the very general symmetrization pro- 
cedure recently proposed by [10]. Since the present method basically exploits the bosonic 
nature of the photons, it can be straightforwardly applied to any particle obeying to the 
Bose statistics. On a more sophisticated technological perspective, the present protocol is 
expected to find applications as a realization of a general, multi-qubit device based on the 
state-symmetrization process to be used in modern Quantum Information networks. [11-13]. 
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Figure Captions: 

Fig. 1. (a) Linear optical scheme for the realization of the general 1 — > M and M — 1 — > M 
Universal Quantum Cloning Machines by a chain of identical symmetrizer beam splitters. 
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(b) Experimental set-up of a 1 — > 3 cloning process. 

Fig. 2. Experimental results of the 1 — > 3 and 2^3 UOQCMs for three input qubits. (a) 
From the upper to the lower row: data corresponding to \{4>4>4>~ L }) an d \{4'4'' L 4'^}) 

components of p^ 3 measured by 3-fold coincidences. From the left to the right column: 
data corresponding to the \H), \H + V), \H + iV) input state p s . (b) From the upper to 
the lower row: data corresponding to \4>4>4>) and \{4>4>4> L }) components of p^ 3 measured 
by 3-fold coincidences. From the left to the right column: data corresponding to the \H), 
\H + V), \H + iV) input state p s . 
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